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(Dated: August 8, 2016) 

The OPERA experiment was designed to search for — > v T oscillations in appearance mode, i.e. 
by detecting the r leptons produced in charged current v T interactions. The experiment took data 
from 2008 to 2012 in the CERN Neutrinos to Gran Sasso beam. The observation of the —> 1 ' 1 - 
appearance, achieved with four candidate events in a subsample of the data, was previously reported. 
In this Letter, a fifth v T candidate event, found in an enlarged data sample, is described. Together 
with a further reduction of the expected background, the candidate events detected so far allow us 
to assess the discovery of iq — > v T oscillations in appearance mode with a significance larger than 5 
a. 


PACS numbers: 


Introduction .— Neutrino flavor transitions due to quan¬ 
tum mechanical mixing between neutrino flavors (v e , iq, 
v T ) and mass eigenstates (iq, zq, v 3 ) were proposed more 
than 50 years ago DU- Several experiments on solar, at¬ 
mospheric, reactor, and accelerator neutrinos have con¬ 
tributed to the understanding of these transitions, re¬ 
ferred to as “neutrino oscillations” Gann. In the at¬ 
mospheric sector, the strong deficit of muon neutrinos 
observed by the Super-Kamiokande experiment in 1998 
was the first compelling observation of neutrino oscil¬ 
lations EH3- This result was later confirmed by the 
K2K jJ5] and MINOS experiments HU- However, for an 
unambiguous confirmation of three-flavor neutrino oscil¬ 
lations in the atmospheric sector, the detection of oscil¬ 
lated neutrinos in appearance mode was required. 

The OPERA experiment has been designed to search 
for Vp —> v T oscillations in appearance mode through 
the detection of the r lepton produced in the v T charged 
current (CC) interactions. It has operated under low 
background conditions and with a signal-to-noise ratio 
as large as about 10. In 2010, a first v T candidate event 
was observed m In 2013, the Super-Kamiokande ex¬ 
periment reported evidence for v T appearance in the at¬ 
mospheric y /t flux with a signal-to-noise ratio of about 
one tenth [l5j. Since 2013, the detection by the OPERA 
experiment of three more candidate events reported in 
Refs. [T3HI5] has allowed us to claim the first observa¬ 
tion of y fl —> y T oscillations in appearance mode with 
a 4.2 a significance [TB]. In 2014, flavor transition with 
high purity in appearance mode has also been observed 
by the T2K experiment in the zq, —» v e channel Ba¬ 
in this Letter, the observation of an additional v T can- 
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didate found in an enlarged data sample is reported. The 
significance of the v T appearance is updated taking into 
account the new observed event and improvements in the 
background evaluation. 

Neutrino beam, detector, and data sample.— The 
OPERA detector at the LNGS underground laboratory 
has been exposed from 2008 to 2012 to the CERN neu¬ 
trinos to Gran Sasso (CNGS) y jL beam [T5|- A total ex¬ 
posure corresponding to 17.97 x 10 19 protons on target 
(POT) resulted in 19 505 neutrino interactions in the 
target fiducial volume. 

The topology of the neutrino interactions is recorded 
in emulsion cloud chamber detectors (ECC bricks) with 
submicrometric spatial resolution. Each brick is a stack 
of 56 1 mm thick lead plates, and 57 nuclear emulsion 
films with a 12.7 x 10.2 cm 2 cross section, a thickness of 
7.5 cm corresponding to about 10 radiation lengths and 
a mass of 8.3 kg. In the bricks, the momenta of charged 
particles are measured by their multiple Coulomb scat¬ 
tering in the lead plates m- a changeable sheet (CS) 
doublet consisting of a pair of emulsion films |20j is at¬ 
tached to the downstream face of each brick. The full 
OPERA target is segmented in about 150 000 bricks ar¬ 
ranged in two identical supermodules (SMs). In each 
SM, the target section is made of 31 walls of ECC bricks. 
Downstream of each target wall, two orthogonal planes 
of electronic target trackers (TTs), made of 2.6 cm wide 
scintillator strips, record the position and deposited en¬ 
ergy of charged particles [21]. A spectrometer, consisting 
of iron core magnets instrumented with resistive plate 
chambers (RPCs) and drift tubes (precision tracker), is 
mounted downstream of each target module. The spec¬ 
trometers are used to identify muons, determine their 
charge, and measure their momentum with an accuracy 
of about 20%. A detailed description of the OPERA de¬ 
tector can be found in Ref. [22]. 

A three-dimensional track in the electronic detector is 
tagged as a muon if the product of its length by the den¬ 
sity along its path is larger than 660 g/cm 2 [23]. An 
event is classified as 1 /r either if it contains at least one 
track tagged as a muon or if the total number of fired TT 
and RPC planes is larger than 19. The complementary 
sample is defined as 0/r. A muon track can be confirmed 
or discarded by measuring its trajectory all along the 
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2008 

2009 

2010 

2011 

2012 Total 

POT (10 iy ) 

1.74 

3.53 

4.09 

4.75 

3.86 

17.97 

Op events 

149 

253 

268 

270 

204 

1144 

lp, events (p M <15 GeV/c) 

542 

1020 

968 

966 

768 

4264 

Total events 

691 

1273 

1236 

1236 

972 

5408 

Detected v T candidates 


1 


1 

3 

5 


TABLE I: Number of events used in this analysis and the 
detected v T candidates for each run year. 


downstream bricks. The momentum-range correlation, 
the energy loss near the stopping point and, eventually, 
the tagging of interaction or decay topologies may con¬ 
tribute to assessing the muonic nature of the track be¬ 
yond the electronic detector performance. 

The analysis described below is extended to all 0/i 
events and to 1 /r events with a muon momentum be¬ 
low 15 GeV/c to reduce the background. The procedure 
starts with the use of the TT hits pattern to select the 
bricks possibly containing the neutrino interaction iM\ . 
These bricks are ordered according to their decreasing 
probability to contain the neutrino interaction vertex. 
The most probable brick (first brick hereafter) is then 
extracted from the target. If the neutrino interaction 
vertex is not found in this brick, it is searched for in the 
next brick in the probability ranking (second brick here¬ 
after). Once the vertex has been located in a brick, a 
surrounding volume of about 2 cm 3 is scanned to detect 
r leptons or other short-lived particle decays [55). The 
details of the event analysis procedure are described in 
Ref. [H]. 

In this Letter, we report the analysis performed on the 
first and second bricks of all of the events recorded by 
OPERA. The event sample is about 15% larger than the 
one reported in Ref. m ■ The numbers of fully analysed 
events are given in Table [I] for each year of data taking. 

The new v T candidate event.— The new v T candidate 
event reported here occurred on August 14, 2012 in the 
second SM, seven brick walls upstream of the spectrome¬ 
ter. As shown in Fig. [l] the activity in the TT is limited 
to the six walls downstream of the vertex brick. The 
event is classified as Op. The visible energy of the event 
is 12 ± 4 GeV. 

A converging pattern of tracks in the CS hints to a 
possible vertex in the brick. Following these tracks in¬ 
side the brick, the neutrino interaction vertex (the pri¬ 
mary vertex) was localised in the 42 nd lead plate from 
the downstream face of the brick. 

The primary vertex consists of the r candidate track, 
which exhibits a kink topology, and a charged particle 
track (PI). The distance of closest approach between 
the r candidate and PI is 0.1 pm, compatible with zero 
within the tracking resolution. In addition to the r lep¬ 
ton and PI, four forward-going and two backward-going 
nuclear fragments pointing to the primary vertex are ob¬ 
served. 


The r candidate decays at a flight length of 
960 ± 30pm into one charged particle which interacts af¬ 
ter crossing 22 plates and can thus be unambiguously 
identified as a hadron. The interaction of the daughter 
particle produces four charged particles and a photon. 
Figure [2] shows the display of the event as reconstructed 
in the brick. 

The difference in angle between the r candidate 
track and the daughter particle track, 0kink, is 90 ± 2 
mrad. The daughter track has an impact parameter of 
83 ± 5 pm with respect to the primary vertex. The z 
coordinate of the decay vertex, Zd ec , measured from the 
downstream face of the lead plate containing the primary 
vertex, is 630 ± 30 pm. A search for nuclear fragments 
has been performed both upstream and downstream of 
the kink vertex up to tan 9 = 3 [25] (with 9 being the an¬ 
gle of the track with respect to the z axis). No fragment 
is found. This result strongly reduces the probability of 
the secondary vertex being due to hadronic interaction. 

The charged particle producing the primary track (PI) 
has a measured momentum of 1.0 ± 0.1 GeV/c. It is iden¬ 
tified as a hadron from its interaction in the downstream 
brick. This, together with the negative search for large 
angle tracks Ea, allows us to rule out the presence of 
a muon at the primary vertex (expected for v^CC re¬ 
lated backgrounds). The linear density of grains along 
the track left by a particle is correlated with the energy 
loss of the particle. The ratio between the grain density of 
track PI and that of the r daughter track is 1.45 ± 0.06, 
to be compared with the 1.38 ± 0.14 expected for a pro¬ 
ton to minimum ionizing particle ratio. Therefore, track 
PI is most likely left by a proton [28] . 

A search for photon conversions possibly pointing to 
the primary and secondary vertices was performed. None 
was found. 

The scalar sum of the momenta of all particles mea¬ 
sured in the brick, p sum , is I 2 I 4 4 GeV/c. The measured 
values of the kinematical parameters and the correspond¬ 
ing predefined selection criteria are summarised in Ta¬ 
ble. n] In the table, p 2ry and p^ v are the momentum 
and the transverse momentum of the decay daughter, re¬ 
spectively, p™ lss is the missing transverse momentum at 
the primary vertex and A<fi T H is the angle between the r 
candidate direction and the hadron direction in the plane 
transverse to the beam direction. The measured values 
of the kinematical parameters of the candidate event sat¬ 
isfy all of the selection criteria for the r —> lh channel. 
The Monte Carlo distributions of the variables and the 
measured values are shown in Fig. [3] 

Signal and background estimation .— The expected 
numbers of signal and background events as well as the 
number of detected v T candidates for each decay chan¬ 
nel are summarised in Table III Assuming Ato 2 3 = 
2.44 x 10 -3 eV 2 [22| and sin 2 20 2 3 = 1, the total expected 
signal is 2.64 ± 0.53 events, whereas the total background 
expectation is 0.25 ± 0.05 events. 
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FIG. 1: Display of the v T candidate event as seen by the electronic detectors in the x-z projection (top panel) and y-z projection 
(bottom panel). The OPERA (right-handed) reference frame is oriented such that the y axis is perpendicular to the hall floor 
and pointing up; the 2 axis is orthogonal to the brick walls and is oriented as the incoming neutrinos. The angle between 
the neutrino direction and the 2 axis projected into the yz plane is 58 mrad. The brick containing the neutrino interaction is 
highlighted in magenta. The solid line shows the direction of the primary track PI (see the text) at its most upstream point 
as reconstructed in the emulsion detectors. 
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FIG. 2: Event display of the fifth v T candidate event in the 
horizontal projection longitudinal to the neutrino direction. 
The primary and secondary vertices are indicated as Vo and 
Vi , respectively. The black stubs represent the track segments 
as measured in the films. 


The numbers of expected signal and background events 
are estimated from the simulated CNGS flux [30]. The 
expected detectable signal events in the Op and 1/i sam¬ 
ples are obtained using the reconstruction efficiencies and 
the v T event rate in the flux normalised to the detected 
interactions. A similar normalisation procedure is 


Parameter 

Measured value 

Selection criteria 

a <Pth n 

151 ± 1 

> 90 

Pt iss (GeV/c) 

0.3 ±0.1 

< 1 

#kink (mrad) 

90 ±2 

> 20 

2 dec (pm) 
p 2ry (GeV/c) 
p^ y (GeV/c) 

630 ± 30 

n±i 4 

l'Oto.4 

[44, 2600] 

> 2 

> 0.6 (no 7 attached) 


TABLE II: Kinematical parameters considered for the r —> lh 
decay channel selection: measured values for the new candi¬ 
date event and predefined cuts are reported in the second and 
third columns, respectively 


also used in the background expectation. The details 
of the signal and background estimation are described in 
Ref. [Til . 

The systematic uncertainty associated with the signal 
takes into account contributions from the limited knowl¬ 
edge of the v T cross section and uncertainties on the sig¬ 
nal detection efficiency. For the signal central value, the 
default implementation for the v T cross-section contained 
in the GENIE v2.6 simulation program is used [31]. A 
10 % model-related systematic uncertainty can be esti¬ 
mated by considering the maximal deviations from the 
central value of the expected number of v T candidates 
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FIG. 3: Monte Carlo distributions of the kinematical variables 
for v T events passing all the location and decay search chain 
with t —> lh decay topology. Red lines show the measured 
values for the candidate event and the corresponding errors. 
Grey areas show the regions excluded by the selection criteria. 


obtained when considering all of the available theoreti¬ 
cal predictions. The only existing measurement of the v T 
cross section is a very low-statistics one by the DONUT 
experiment [52]. Owing to the fact that the v T signal ex¬ 
pectation is calculated by using location efficiencies de¬ 
termined from the lfj, and 0 /i data samples, this value 
is at first order insensitive to systematic effects on ef¬ 
ficiencies up to the primary vertex location level. Fur¬ 
ther confidence on the global efficiency estimation is ob¬ 
tained by considering the charm data sample for which 
good agreement is found between the 50 observed events 
and the expectation (54 ± 4) provided by the neutrino- 
induced charm production cross section and the detec¬ 
tor simulation j25J ■ Additional uncertainties on the 
number of expected v T candidates arise from the experi¬ 
mental knowledge of 6*23 and Ato | 3 (10 %), and from the 
uncertainty in the efficiency for tagging r lepton decays 
(15%). The latter contribution arises from the statistical 
error of the sample of iA c events with charm production 
which was used for validation. The CNGS flux uncer¬ 
tainty plays a minor role since the expected number of 
v T events is determined from the detected interac¬ 
tions used as a normalisation sample. The simulation 
of the kinematical properties of the final state was per¬ 
formed using the NEGN generator [51] , which takes the 


polarisation of r leptons into account (r decay library 
TAUOLA [35]). The associated systematic uncertainty 
on the expected number of r decays in all channels is 
estimated at the level of a few percent [3B]. The total 
systematic uncertainty on the expected signal is then set 
to 20 %. 

The main processes contributing to the background 
for the v T appearance search are charmed particle de¬ 
cays, hadronic interactions and large-angle muon scat¬ 
tering (LAS). The corresponding contributions are es¬ 
timated by simulation studies validated with real data 
samples. Using the measured sample of CNGS r^CC in¬ 
teractions with charm production, the uncertainty on the 
charm background has been estimated to about 20% [25]. 
This includes a contribution from the experimental un¬ 
certainty on the charm cross section ( 8 % [32]), the hadro- 
nisation fraction ( 10 %), and the statistical error of the 
CNGS charm control sample (15%). Hadronic back¬ 
ground has an estimated uncertainty of 30% from data- 
driven measurements of test-beam pion interactions in 
the OPERA bricks |37| . 

With respect to what was reported in Ref. El, 
an additional improvement in the estimation of the 
LAS background in the r —>■ /i decay channel has 
been achieved (381- The LAS rate is estimated us¬ 
ing a GEANT4-based simulation implementing a mixed- 
approach algorithm with ad hoc modifications to take 
into account the effect of the nuclear form factor at the in¬ 
volved transferred momenta (of the order of a few fm -1 ). 
A Saxon-Woods charge density is assumed with param¬ 
eters derived from fits to data. Scattering off individual 
protons is also taken into account. The simulation is 
benchmarked on experimental data including scattering 
of 2 GeV/c muons on a 12.6 mm lead target, 7.3 GeV/c 
and 11.7 GeV/c muons on a 14.4 mm thick copper target 
and 0.512 GeV/c electrons on a 0.217 mm lead target [391 
HD- From this study, it follows that the number of LAS 
background events that satisfy the r —> [i selection crite¬ 
ria amounts to [ 1.2 ± 0 . 1 (stat.) ± 0 . 6 (sys.)] x 10 ~ 7 /z^ c 
interactions, well below the conservative value considered 
in our past publications. 

Results .— In this analysis, the observed number of 1 ',- 
candidates rii for each individual r decay channel i is 
considered as an independent Poisson process with ex¬ 
pectation fiSi + bj . The expected signal and background 
events, Si and bi respectively, are taken from Table |III| 
the signal strength factor /j, is a continuous multiplicative 
parameter for the expected signal. The background-only 
hypothesis corresponds to /r = 0 , and the nominal signal 
to /i = 1 . 

The significance of the observed v T candidates is eval¬ 
uated as the probability that the background can pro¬ 
duce a fluctuation greater than or equal to the observed 
data. Two test statistics are used for the computation; 
in both cases, the test statistics values of the observed 
data are compared with sampling distributions obtained 
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Channel 

Expected background 

Expected signal 

Observed 

Charm 

Had. reinterac. Large p scat. 

Total 

r —> lh 

0.017 ±0.003 

0.022 ± 0.006 

0.04 ±0.01 

0.52 ±0.10 

3 

t — > 3h 

0.17 ±0.03 

0.003 ±0.001 

0.17 ±0.03 

0.73 ±0.14 

1 

t —y p 

0.004 ±0.001 

0.0002 ± 0.0001 

0.004 ±0.001 

0.61 ± 0.12 

1 

t — » e 

0.03 ±0.01 


0.03 ±0.01 

0.78 ±0.16 

0 

Total 

0.22 ±0.04 

0.02 ± 0.01 0.0002 ± 0.0001 

0.25 ±0.05 

2.64 ±0.53 

5 


TABLE III: Expected signal and background events for the analysed data sample. 


with pseudoexperiments. 

The first test statistics is based on the Fisher’s method. 
For the background-only hypothesis (i.e., p = 0), the p 
values pi of each individual channel (calculated as the 
integral of the Poisson distribution for values larger or 
equal to the observed number of candidates) are com¬ 
bined into an estimator p* — n, Pi [121 03]. By com¬ 
paring the observed p(] ata with the sampling distribution 
of p*, a (one-side) significance of 5.1 standard deviations 
is obtained, corresponding to a background fluctuation 
probability of 1.1 x 10 '. 

The second test statistics is based on the one-sided 
profile likelihood ratio X(p) [29] . This test statistic is 
used to quantify the discrepancy between the data and 
a certain hypothesised value of /i. The significance, the 
level of disagreement between the observed data and the 
p = 0 hypothesis, is computed by comparing Ad a ta(/^ = 
0) with the corresponding sampling distribution of A(/i = 
0). The likelihood, which includes Gaussian terms to 
account for the background uncertainties, is 


The 90% confidence interval for Am | 3 has been es¬ 
timated with three different approaches using the pro¬ 
file likelihood ratio, the Feldman-Cousins method, and 
Bayesian statistics. Assuming full mixing, the best fit 
is A777-23 = 3-3 x 10 -3 eV 2 with a 90% C.L. interval of 
[2.0, 5.0] x 10 -3 eV 2 , the differences among the three 
methods being negligible. 


Conclusions This Letter reports the analysis of a 
data sample including the first and the second most prob¬ 
able bricks for all runs, with a corresponding increase of 
the statistics of about 15% with respect to Ref. [16]. In 
this enlarged data sample, a fifth r neutrino candidate 
has been found. Furthermore, a revision of the back¬ 
ground estimate in the muonic decay channel has been 
performed. Given the low background level and the ob¬ 
served number of v T candidate events, we report the dis¬ 
covery of a v T appearance in the CNGS neutrino beam 
with a significance of 5.1 a. 


C = [QPoisson(? 7 i |/is i +/3 i )Gauss(/3 i |6 i ,CT i , i ), (1) 


i=1 


where cp,. is the background uncertainty for channel i 


(from Table III) and fit are the background parameters 


Gaussian modelled. Two different implementations of 
the method, one based on a custom code and the other 
one based on RooStats [44] . have been used with both 
giving a significance of 5.1 standard deviations. 

A simple compatibility test of the observed data with 
the expectations from the neutrino oscillation hypothesis 
(p = 1) is given by the best-fit signal strength at 90% 
C.L., p = 1.8^[' 3 , which is consistent with unity. An¬ 
other test was made by performing pseudoexperiments 
to sample the distribution of the data assuming p = 1 
and taking into account the uncertainties on the expected 
signal and background. The probability of data being 
less likely or equal to the observed ones is 6.4%. If we 
consider the total number of v T candidates regardless of 
the distribution into decay channels, the probability of 
observing five or more candidates with an expectation 
of 2.64 signal plus 0.25 background events is 17% from 
Poisson statistics. 
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